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Abstract 

The effective potential of scalar quantum electrodynamics with N flavors of 
complex scalar fields is studied, by performing a self consistent 1/N expansion 
up to next to leading order in Starting from the broken phase at zero 

temperature, the theory exhibits a phase transition to the symmetric phase at 
some finite temperature Tc. We work in general covariant gauges and demonstrate 
the gauge invariance of both, the critical temperature Tc and the minimization 
condition at any finite temperature T. Furthermore, the only minimum of the 
potential is at zero scalar vacuum expectation value for any temperature T > Tc 
and varies continuously to nonzero values for temperatures below Tc, implying the 
existence of a second order phase transition. 
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1 Introduction 



In the last years, there has been a great interest in understanding the nature of the finite 
temperature symmetry restoration phase transition in the Standard Model (SM) ^ - 
|T2| . Recently, the main motivation arose from the realization that the rate for baryon 
number violation in the SM at finite temperature is much larger than what it was orig- 
inally assumed - In fact, once the temperature is much higher than the height of 
the energy barrier between the different topologically distinct vacua the exponential 
suppression of the anomalous baryon number violation processes, present at low temper- 
ature, disappears 0-^]. Thus, at sufficiently high temperatures, the anomalous baryon 
number violation processes could wash out the primordial baryogenesis generated, in 
some grand unified scenarios, at the very early stages of the universe |]T4[. In such cases, 
the generation of baryon number could have occurred through non-equilibrium processes 
in the transition from the symmetric to the broken SU{2)l x f/(l)y phase 0. Quite 



generally, models of electroweak baryogenesis [0 fulfill Sakharov's three requirements 
|16| for the generation of baryon asymmetry. Besides the nonconservation of B the un- 



derlying dynamics must involve nonequilibrium processes and violate CP. The violation 
of CP is present in the SM and its extensions, although in the SM it is probably too weak 
in order to generate the required baryon asymmetry [T^ . To satisfy the nonequilibrium 
condition, these models rely on the presence of a first order phase transition, sufficiently 
strong in order to imply the suppression of baryon number violation processes in the 
broken phase. 

The nature of the electroweak phase transition is a fundamental question which 
needs to be carefully explored and requires a detailed study of the effective potential of 
gauged scalar theories at finite temperature. A one loop analysis of the abelian theory 
was first done by Dolan and Jackiw 0. The one loop effective potential shows the 
existence of a first order phase transition, whose strength depends on the values of the 
electroweak gauge and the scalar quartic coupling constants However, the naive 
one loop definition of the effective potential is known to have serious infrared problems, 
such as the appearance of a complex effective potential even at temperatures above 
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the critical temperature, Tc, at which the transition from the broken to the symmetric 
phase takes place. Moreover, it is found that the critical temperature itself is a gauge 
dependent quantity. It was realized in Ref.0 that both - the infrared and the gauge 
dependence - problems do not appear, if only the dominant terms in the high temperature 
expansion are kept. Within such an approximation, a definition of a gauge invariant 
critical temperature can be achieved by assuming the existence of a second order phase 
transition. However, in order to study the nature of the phase transition, higher order 
terms in the effective potential need to be computed. 

Quite recently, an improved analysis of the one loop effective potential has been 
considered by several authors -[|12|. The improved one loop effective potential is 
generically obtained by the inclusion of the finite temperature one loop corrections, at 
zero external momentum, to the boson propagators. When this process is carefully done, 
in order to avoid double counting of diagrams, it is equivalent to the inclusion of the so 
called ring diagrams in the one loop computation |T^. Such resummation of diagrams 
solves several problems inherent to the one loop expansion and leads to a first order 
phase transition, although weaker than that one found within the one loop expansion 
||1 1|| . However, in the pure scalar case the ring expansion breaks down at temperatures 



close to the critical one, |T — ^ ATc, with A the scalar quartic coupling Hence, it 
can not be used to analyze the behavior of the effective potential at T ~ Tc. In addition, 
after the inclusion of the gauge fields, for values of the squared gauge coupling of the 
order of the quartic coupling it is not possible to conclude that the phase transition is 
first order based solely on the ring expansion. This is due to the fact that, within the 
ring approximation, and for ~ A, the difference between the transition temperature 
To, at which the curvature at the origin vanishes, and the critical temperature Tc is 
of order g'^T^, but the ring expansion breaks down for |T — T^ < g^T^ [Q. Thus, an 
analysis beyond the improved one loop approximation is necessary in order to study 
the nature of the phase transition for g'^/X — 0{1), which is the phenomenologically 
interesting region in which the Higgs mass is of the order of the gauge boson mass. 

As it has been explained in the seminal work by Dolan and Jackiw [0, for the pure 
scalar case a self consistent resummation of the diagrams contributing to the effective 
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potential can be done, by extending the theory to include N flavors of self interacting 
scalar fields and performing a 1/N study of this model. The 1/N expansion provides 
information which does not appear at any fixed order of perturbation theory and, hence, 
it is a more appropriate tool to study the nature of the phase transition, which is in 
itself a nonperturbative phenomenon. Studying the behaviour of the effective squared 
mass, by doing a self consistent 1/N expansion, Dolan and Jackiw showed that the 
infrared problems inherent to the pure 0(N) scalar theory disappear. The explicit form 
of the effective potential at finite temperature was, however, not derived in this analysis, 
and it was studied in Refs. and [|T^. A similar resummation of diagrams has been 
recently performed in Ref . |jl9| , while an analysis of the component - 0^ theory, by a 
method based on average fields, has been considered in Ref. . 

The validity of the finite temperature computations in the scalar 0(N) linear model 
was questioned, when it was realized that, in the continuum limit, the effective scalar 
potential at zero temperature was either complex, unstable or did not allow the sponta- 
neous breakdown of the 0(N) symmetry ||21|] - This issue was clarified by Bardeen 



and Moshe Moshe ||2J] , who explained that the problems appearing in the large N com- 
putations were not associated with the 1/N expansion but with a more fundamental 
nonperturbative property of the 0(N) scalar model in four dimensions, which is the 
issue of triviality. In fact, for positive bare quartic coupling, the continuum limit of 
the four dimensional 0(N) theory is trivial, in the sense that the renormalized coupling 
goes to zero when the cutoff is removed. This implies that, to work with a nontrivial 
theory, one has either to take a negative bare quartic coupling, rendering the theory 



unstable at finite temperature |]2^, or treat the theory as an effective one, by keeping a 
finite effective cutoff. For a large effective cutoff, the 1/N expansion at zero and finite 
temperature can be defined in a consistent way, avoiding the problems mentioned above. 
Thus, in the following, we shall assume that we are dealing with an effective theory valid 
up to an energy scale of the order of an effective cutoff A, at which, quite generally, new 
physics should appear. 

It must be mentioned, that the nature of the phase transition for scalar QED with A^ 
fiavors and a scalar Higgs heavier than the gauge boson, have been also studied by using 
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a 4 — e expansion |25|-||2^. Within such approximation, the theory may be analyzed by 
performing a renormahzation group study for small e. The high temperature theory may 
be treated as an effective three dimensional theory and thus, to obtain physical results, 
the expansion parameter must be taken to be e = 1. In this framework, the critical value 
for the number of flavors, Nc, below which the phase transition becomes first order, was 
found to be surprisingly large, N^. = 183. Based also on this analysis, Halperin, Lubensky 
and Ma |^ predicted a first order phase transition for the smectic-A to nematic phase 
transition in liquid crystals. However, when contrasted with experiment, this transition 
was found to be second order [E^]. In addition, for e 1, a (2 + e) expansion gives 



results which are at variance with those ones coming from the (4 — e) expansion |25||,P8 
In particular, in the case in which the Higgs particle is heavier than the gauge boson, 
the (2 + e) expansion indicates a much lower critical value of A^, if any, for the three 
dimensional theory ||28 



This paper is organized as follows: In section 2 we present the model, a scalar theory 
with N flavors of complex scalar fields coupled to an abelian gauge field. In section 3 
we derive the general expression of the effective potential up to next to leading order 
in 1/N. In section 4 we concentrate on the leading order results, and we show that a 
second order phase transition takes place in this case. In section 5 we carry out the 
analysis of the critical temperature and the nature of the phase transition up to next 
to leading order in 1/N. We start discussing three related aspects of the theory, which 
are the renormahzation of couplings, the triviality of the theory and the absence of 
tachyons. Then, we analyze the radiative corrections at finite temperature to finally 
perform the minimization of the effective potential up to next to leading order in 1/A^. 
We demonstrate that a gauge invariant expression for the position of the minimum is 
obtained. In addition, the position of the minimum is not modified by next to leading 
order corrections and, hence, the phase transition remains second order. A comparison of 
our large A^ results with those previously derived within the improved one loop method 
is performed in section 6, where we analyze the source of the discrepancy between the 
results obtained in both approaches. We reserve section 7 for our conclusions. 
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2 The Model 



The Lagrangian density for a scalar theory with N flavors of complex scalar fields in 
interaction with an abelian gauge field is given by 

In the above expression, a is a flavor index taking values from 1 to A^, 0" is a complex 
scalar field, |0p = (pi(pa, is an abelian gauge field, and is the covariant derivative 
associated with it. In addition, we have included a gauge fixing term with a being the 
gauge fixing parameter (a = is the Landau gauge). Note that, at vanishing gauge 
coupling, this model reduces to the 0(2N) linear sigma model. 

Following the method first proposed by Jackiw , to compute the effective potential 
V{(f)) for this theory, we shift the fields 0a — ^ 4>a + (pa, where (pa is the vacuum expectation 
value of the scalar field and (pa stands for its quantum fluctuations. The Lagrangian must 
be truncated by dropping all terms which are linear in the (pa quantum fluctuations. The 
new Lagrangian, C{(pa + (pa,Afj), may be decomposed into a term quadratic in the 
and (pa fields and an interaction term. The effective potential is given by the sum of the 
tree level potential plus the one loop contributions, obtained from the bilinear part of 
the shifted Lagrangian, plus higher order loop contributions, given by the sum over all 
connected one-particle irreducible vacuum graphs for the theory described by the shifted 
Lagrangian density C{(j)a + (pa, A^). 

The quartic interactions pose a disadvantage for the diagrammatical analysis of the 
theory. Therefore, following Refs. |2T| - we shall eliminate the scalar quartic term 
through the introduction of an auxiliary field x- The modified Lagrangian reads 



3 . r 

C{(Pa,A^,X)= ^(l)a,A^) + -^ 

LA 



^"6 



A /, , ,o 2 



(2) 



The Euler-Lagrange equation for x is a constraint equation relating x with the scalar 
field 0. In fact, the effective potential for the modified theory V"(x, 0a) reduces to the 
one of the original theory if x satisfies the requirement 



which defines the vacuum expectation value x as a function of 0. 

A self consistent expansion can be defined if the gauge and quartic coupling 
constants g and A, respectively, depend on so that in the limit of ^ oo, g'^N and 
AA^ go to constant values. Hence, for the purpose of our study it is better to define 
the new coupling constants e = gy/N and A = AA^. Once the field x is introduced the 
Lagrangian density of the theory reads 



+ ^A,A^<Pa,<Pa, + ^-^x' -liW - V') , (4) 



2N~ ^ ~ """""" ' 2A " 2 

where i = 1,2, (f)a,i and 0^,2 are the real and imaginary parts of the field (pa, respectively, 
and summation over a and i is understood. Although the introduction of the field x does 
not alter the dynamics of the full theory, it does lead to a new perturbation series, in 
which the 1/N expansion has a simpler diagrammatic interpretation. After considering 
the shifted fields, which now include the shift for the auxiliary field x, 



'N<Pa + 0a, X^X + X (5) 

and ignoring the linear terms in the quantum fiuctuations x and 0a, the shifted La- 



- ^ {d^A^ f + i-A^A^^a,k, - eA^^e'^4>a,,d^(f)a, 
la I 

(l)a,jd^(l)a,i + —A^'Af,(f)a,i(f)a,t + ^= A^' A^Cpa.iCpa,! + -'9/.0a,»<9''0a,i 

/A V A" 2 

+ - ^X0a,i0a,i - ViV0a,iX0a,i - ^0a,i0a,iX - '^iree, (6) 

where the tree level potential is given by, 

V^ree = -^-§^x' + ^x{\^?-v'). (7) 



rang 


ian reads 


C 








e 








Vn 




3N 



In the above, we have followed Root [g^ in rescaling the vacuum expectation value 
0a and, consistently, f by a factor y/N. This simplifies the 1/A^ counting. In fact, if 
this rescaling were not done, the vacuum expectation value would be of order N in the 
natural scales of the theory, 0^ = A^F(/i^,T^), where F is a function of the bare mass 



/i^ = — the temperature T and the couphngs of the theory. This extra factor of 
N would have to be properly considered while doing the expansion at next to leading 
order. The physical results, of course, would remain unchanged. 

The above Lagrangian may be decomposed into a quadratic and an interacting part. 
The expression of the quadratic Lagrangian reads 

1 



J — -Ai^ 

■'quad 2 



OAT 

+ -^X' - VN4>a,iX<Pa,i - eA^e'^,,^d,(t>a,i (8) 

from which it is easy to deduce the inverse propagator matrices relevant for the com- 
putation of the effective potential. The interacting part of the Lagrangian £/, on the 
other hand, determines the vertices of the boson fields interactions. For simplicity, we 
will assume in the following that 

01,1 = k (9) 

while all other vacuum expectation values are zero. In fact, it is not important which 
state is chosen as the ground state of the theory. Independent of such a choice the 
original global U(N) symmetry is spontaneously broken to U(N-l) and 2N-1 Goldstone 
bosons are generated, one of which is eaten by the gauge field through the usual Higgs 
mechanism. 

We obtain two decoupled inverse propagator matrices. One that mixes the gauge 
field An with 0i_2, 0i,2), and the other one which mixes the auxiliary field x with 

01,1) -^(x^^i.i)- Thus, in momentum space representation, the nonzero components of 
the inverse propagator matrices read. 



SabSijik'^ - x) 

3N 

X 



(10) 
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Considering the above expressions it is straightforward to compute the determinant of 
the inverse propagator mass matrices, which in Euchdean space are given by 



det7W(A 



(P + x)(A;2 + e202)3 



a 



detM{x, 



-3N 

IT 



X 



(12) 



Furthermore, the propagators of the scalars and gauge bosons, which are relevant for 
the computation of the effective potential up to next to leading order in may be 
then obtained from Eqs. (p!OD , ( p!TD and ([T^), and in Euclidean space are given by 

1 k^,k^ a{k'^ + x) 



iD 



fJ,U 



5 



fJ,U 



k^ k^ + e 



2A2 



A;2 



fc^ + xk"^ + ctxe^cf)'^ 



iD 



(k' + x) \/3N 



XX 



F + X + A02/3 



(13) 
(14) 



3 The Effective Potential 



As we mentioned above, in order to analyze which Feynman diagrams contribute to 
the effective potential up to next to leading order in 1/A^, it is necessary to consider the 
behaviour in 1/A^ of the propagators and of the interaction vertices. From Eqs. ([T0D-(p!4D 
we observe that a factor 1/A^ is associated with D^^ and a factor 1/\/N is associated 
with -Di,i;x- Moreover, a factor of order is associated with any closed loop of (pa fields, 
due to the summation over all possible internal fields. From the interaction part of the 
Lagrangian in Eq.(|^), we observe that a factor is associated with the derivative 

coupling of the scalar fields (pa to the gauge bosons, while a factor 1/A^ appears in any 
quartic vertex A'^(j)a,i(f>a,i- Note that the vertex X(pa,i<Pa,i has, instead, a factor 1 in the 
coupling. 

The only diagram contributing to the potential at leading order in 1/A^ 1113 'EH' ^i-o., 
is the closed loop involving the 2(N-1) Goldstone bosons which do not mix with either 

or X, 

K„, = H<^jfl„(fc^ + x). (15) 
9 



We perform the finite temperature computations in the imaginary time formahsm: After 
a Wick rotation to Euchdean space we impose periodic boundary conditions on the time 
direction of length L = = fi. For simphcity of notation, we have defined 



with ujn = 27rn//3. 

The computation of the effective potential at next to leading order in involves 
two type of contributions. There are one loop contributions, involving the two nontrivial 
propagator matrices. 



and there are those coming from multiloop one particle irreducible graphs, depicted 
in Fig.l. Observe that any multiloop contribution involving the mixing of x ^^id A^, 
like those depicted in Fig. 2, vanish in the regularized theory. After some work, the 
resummation of all the multiloop diagrams contributing to next to leading order in 
gives 




(16) 




(17) 




(18) 



which can be rewritten as 



n. 



n—loop 
i.l.o 




(19) 



where 




(20) 



and 




(21) 
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is the finite temperature vacuum polarization contribution. Furthermore, using the 
properties 



detM{(j)i,2, A^) = [k^ + x) det [-iD^^l 



(22) 



and 



detMi^,,i,x) = {k' + x)D-i, (23) 
we arrive to the final formal expression for the effective potential up to next to leading 



order in 



Vtree + Vl,o. + V^j[^^ + .o""^ 



X + yX (I 



2A 



^In {k^ + x 



2 Jk 



In det 



iDj{k, X, (f)) - e^Yi^u{,k, x) 



+ - / In 

2Jk 



{k^ + x){i + \B{x,k^)/z) + M^/?, 



(24) 



A;2 + X 

As we discussed above, the expression for the effective potential of the original theory 
may be obtained from the above equation by imposing the condition 

dV{^,x) 



dx 



0, 



(25) 



which determines the expression of x as a function of cfP'. Furthermore, as it was first 



noticed in Ref . |^ , in order to obtain the effective potential up to next to leading order 
in it is sufficient to solve the gap equation, Eq. (P^D, up to leading order in 
In fact, calling x(0^) the exact solution to the gap equation and x{4>^) the leading order 
solution, then. 



X(0^) = 0{l/N). 



Expanding the effective potential around x(0^) = x{4>^)^ we have 

dV{^,X) 



V(0,X) = V^(0,X) + 



dx 



(X-X) + 0{IIN). 



(26) 



(27) 



However, since x(0) is the solution to the gap equation at leading order, it follows that 
<9VJ.o. (0, x) /f^X = 0- Therefore, the second term in the above equation is also of order 
1/N, 

\/(0, m) = V{^, xm + 0{1/N). (28) 
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Thus, in order to compute the effective potential up to next to leading order in we 
just need to consider Eq.(p4D with x given by its leading order expression, = 

Moreover, for the extent of this work we shall concentrate in computing the location 
of the extrema of the effective potential, which are derived from the condition, 

dV{4>,x) 



0. 



(29) 



In fact, the number and location of the extrema of the effective potential provide suf- 
ficient information to study the nature of the phase transition. In addition, as we will 
show below, a gauge independent expression for the solutions to Eg . (|29|) is found. Ob- 
serve that if Eq.(^) is not fulfilled for any nontrivial value of the scalar field, the only 
minimum would be at the origin and the global U{N) symmetry, together with the local 
U{1) symmetry of the theory would be preserved. From Eg . (plD , we find the relation 
dVio./dcf)'^ = Nx/2, and recalling the fact that dVio./dx = 0, we obtain 



dV 



iVx(0 '^ 
2 

/ d 



+ 



+ 



dx d 



In det 



+ - / In 

2Jk 



2 ■ d(t)'^dx) 12 J k 

(fc2 + x)(l + Ai?(x,fc2)/3) + A0V3 



l{k, x,0) - e^Tl^,^{k, x) 



(30) 



k"^ + X 

which is the formal expression from which, after proper integration and renormalization 
procedures, we shall be able to determine the critical temperature and the order of the 
phase transition up to next to leading order in 



4 Critical Temperature and Nature of the Phase 
Transition at Leading Order in 

In section 3 we derived the expression of the effective potential up to leading order in 
to be. 



v{xA) 



2X. 



-2 



N 
~2 



+ —X (<; 



N l^ln (A;2 + x) 



(31) 



where we have introduced the subscript zero to denote unrenormalized quantities. Thus, 
at leading order in 1/N, the gap equation, which determines the expression of x as a 
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function of dP', reads 



dV_ 

dx 



3N _ N 



0. 



(32) 



Ik k"^ + X 

In order to explicitly evaluate the temperature dependent integrals, we shall use the 
relation |^ 



f{k) + / ^ 



°°+*^ dk^ 

oo+ie 



d^k 



f{ko,k) + f{-ko,k) 

{exp{-iPko) - 1) 



(33) 



(27r)4'' ' ' J-oc+ie 27r J (2%)^ 

which allows a decomposition in the zero temperature and the finite temperature con- 
tributions - the first and second terms in the right hand side of the above equation 
respectively. For the particular case of the integral appearing in Eq.(^2|), we obtain, 
1 r d^k 1 f d^k 1 



+ 



kk' + X J (2vr)4 (P + x) J (2vr)3 ^ 



+ X 



exp (3Jk^ + x]-l 



(34) 



where the fco dependence in the last integral was evaluated by performing a contour 
integration in the complex /cq plane. At large temperatures, T 3> x, the above integral 
may be expanded in a way first derived by Dolan and Jackiw 



d^k 



1 



Vx 



k"^ + X 



exp Pyk"^ + X 



+ 



12/?2 4vr/? 167r2 



X 



+ h.o.{P'x) 
(35) 



where c = 167r2exp(l — 27), with 7 ^ 0.577. 

Moreover, the zero temperature contribution is quadratically divergent and needs to 
be regularized. Introducing a momentum cutoff A, so that 



d^k 



A' -xln 



'A^ 



(36) 



{27r)^k^ + x) 167r2 L ^ \x 
we then absorb the quadratic and logarithmic divergences by defining the renormalized 
quantities 



and 



1 
A 



A^ 



1 1 



In 



A^ 

M2 



(37) 
(38) 
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where is a renormalization scale. Observe that, for any positive value of the unrenor- 
malized quartic coupling Aq, the renormalized coupling A vanishes at A — oo, and the 
theory becomes trivial in the continuum limit. If the unrenormalized quartic coupling is 



taken to be negative, the theory becomes unstable [^. Therefore, the only consistent 
definition of the theory is when it is considered as an effective theory with an effective 
finite cutoff A. 



We can further define a temperature dependent coupling by 



1 _ 1 1 
A^ ~ A ~ 48^ 

to finally rewrite the gap equation as 

At 12' ^ 



In 



cT2 
M2 



(39) 



X 



6 



6p^ 



(40) 



The above expression gives a quadratic equation in y/x, which can be easily solved to 
find 0, p 



X 



, ( \ +^ 

2471(3 \ \24n(3j 6 



(41) 



where the plus sign has been chosen in order to obtain positive values of y/x, which is 
a precondition for the validity of the gap equation derived above. 

As we have shown in section 3, the minimization condition up to leading order in 
1/A^ reads, 

0. (42) 



2^ . Nxii') 



Considering the explicit expression for x(0^)) given in Eq. (PT|) , this yields the relation 

1 



6/32 



0, 



which determines the value of at the minimum. Eq. 
atures above the critical value 

^ 6v\ 



rp2 — ~ 



(43) 

has no solution for temper- 

(44) 
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and, hence, the minimum is at 0mm = in such region of T.Q On the other hand, at 
temperatures below the critical one, the vacuum expectation value is given by 

1/1 1 \ 



(45) 



Equation (|^) shows that the vacuum expectation value varies continuously from zero 
to nonzero values, signaling the presence of a second order phase transition within the 
leading order in 1/A^ solution of the model. It is worth to remark that, for values of 
0^ at the left of the minimum, 0^ < {v"^ — 1/6/3^), -y/x takes negative or even complex 
values, showing that, as was first discussed by Coleman, Jackiw and Pohtzer the 
effective potential can not be defined in a sensible way in that region. 

5 Critical Temperature and Phase Transition: Anal- 
ysis up to Next to Leading Order in 1/N. 

In terms of the renormalized quantities defined in section 4, and within the framework 
of the high temperature expansion, we can rewrite the effective potential up to next to 
leading order in 1/A^, Eq. (p^ , as 



T/ri 3iV _ iV _ 2\ ^ Nx Nx'^^ 

y{<l>,X) = -TT^X -V + 



+ 



\ Ik ^^^^^ H^Mi'^^^' ^' ~ e^n^^(/c, x) 

{e + X) {I + \oB{x,k)/3) + Xo4>VS 



+ - fin 

2Jk 



(46) 



k'^ + X 

As we shall explicitly show in section 5.1, B{x,k^) and Il^jy (x, fc^) are logarithmically 
divergent quantities. However, these logarithmic divergences can be naturally absorbed 
in the renormalization of the gauge and quartic couplings entering in the integrands of 
Eq.(^). The remaining divergences, arising from the zero temperature contributions to 
the integrals, may be absorbed in the renormalization of A and f ^ at next to leading 
order in 1/A^. 

^ Observe that there is a difference in a factor 2 with respect to the results of Refs. due 
to the presence of 2N, instead of N, real scalar bosons in the theory. 
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From Eq.(30) the minimization condition is given by 
1 



X 



N 



d dy d 



62 dx 



In det 



-iD l{k, X, <j)) - eon^^(/c, x) 



In 



{e + x) (l + Aoi?(x,A;)/3) + Ao0V3 



(47) 



x=o 



k'^ + X 

Thus, since the minimum is obtained for x — ^{^1^)-, we can safely set x = in the 
right hand side of Eg. (|^ , under the assumption that no infrared divergences arise in 
this process. 



5.1 Tachyon Poles and Triviality of the Gauged 0{2N) Model 
at Zero T 

At zero temperature, it is possible to compute the radiative corrections B{p,x) and 
^nu{p, x) by making use of a gauge invariant regularization scheme, like Pauli - Villars. 
This gives the result 

Bt=o{p\x{4>)) 



- 2 



while 

^^=o(p^x(0')) 



'4x+_p' 



X I 

2\ 1/2 



In 



X U 4p2 ] +2 



+ 2 



487r2 
'^X±^ 



, A2 2 
In 

X 3 



(4^ 



(49) 



+ p' 
4p2 



2\ 1/2 



In 



^ H + 

x\\ 



2\ 1/2 



The divergence associated with B{p'^,x) may be absorbed in the renormalization of A 



22| . In fact, the expression 

F{x,p') ^ (p' + x) (3/Ao + Br=oix,p')) + ^\ 



appearing in the potential, Eq.(0), is proportional to the inverse propagator of the 
massive scalar particle a = and may be rewritten as 

Fix.p") = {p' + x) (3/A + Bt=o{x,p')) + 0', (51) 
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where 



5t=o(p,x(0')) 



167r2 



In 



M2 



4p2 



X J 

2\ 1/2 



In 




4x + P^ 
4p2 



1/2 



+ 1 



1 



(52) 



and A is the renormahzed quartic couphng at leading order in which is given in 

Eq.(p 



Observe that, as first noticed in Ref. |21], the renormahzed theory seems to be 
spoiled by the presence of tachyons. In fact, since we are working in Euclidean space, a 
tachyon pole in the 0i^i two point function will appear if, for some positive value of p'^, 
PiXiP^) = 0- Considering the limit x = = 0, this implies an equation for p^, which 
reads 

?s 1 / A/r2 \ 1 

0. (53) 



3 
A 



For a finite value of the renormahzed coupling A, the above equation has a solution at 
sufficiently large values of p^, which could be interpreted as the presence of a tachyon 
in the spectrum. However, rewriting Eg . (|53|) in terms of the unrenormalized quartic 
coupling, we obtain, 

1 /A2\ 1 

0. (54) 



3 1 
A^^"^ 16^ 



Therefore, for any Aq > 0, the potentially dangerous pole appears at momentum above 
the cutoff scale [ |24|| , 

^2 A 2 f /i^^ 



p2 ~ a2 exp (487rVAo) • (55) 

Hence, it has no physical consequences. (Observe that at nonvanishing values of and 
X the pole would be at even larger values of p^.) The apparent discrepancy between the 
results obtained in terms of the renormahzed and unrenormalized couplings is related 
to the issue of triviality. In fact, while naively analyzing the existence of a tachyon pole 
using Eq.(l^), one overlooks the fact that the renormahzed quartic coupling A goes to 
zero in the continuum limit A — oo, and for this reason, a finite effective cutoff is needed 
in order to have a finite renormahzed coupling. 

An analogous situation occurs in the gauge sector of the theory. In fact, calling 



iD 



p^p" 

p2 



D 



■' 9 



(56) 
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and 

nr=o(x,p') = [s'^'' - ^) ^T=oix,p'), (57) 
a tachyon pole in the current - current correlation function would appear if 

G [p\ x) = \Dr' (x, 0, - nT=o (x, (58) 

had a zero in Euclidean space. At x = = 0, the condition (^(p^, x) = leads to the 
equation 



1 



0. (59) 



el 48772 

Apart from a change from the gauge to the quartic coupling, the above equation is 
equivalent to that one found in the previous case, Eq. (|5^) . The would be tachyon pole 
appears at momentum above the physical cutoff of the theory 

~ A^exp(487rVeo) (60) 

and, hence, has no physical implications. 

The renormalization of the gauge coupling proceeds in similar way to that one of the 
quartic coupling. The logarithmic divergences in IIj.'^q are absorbed in the definition of 
the renormalized gauge coupling e, which is given by 

1 1 1 , /A'\ 

^ = ^ + 48^^H^J- ^ ^ 

We can therefore define the renormalized vacuum polarization, 11^=0; by replacing 
ln(A2/;^) by \n{M'^/x) in IIt=o- Moreover, the expression of the renormalized gauge 
coupling, Eq.(^), shows the triviality of the gauged theory in the continuum limit. 

5.2 Radiative Corrections at Finite Temperature 

We shall now analyze the structure of B{p, x) and Il'^'^{p, x) at finite temperature, which, 
by making use of Eq. (|33|) , may be decomposed into zero temperature contributions, Bt=o 
and HtLo^ already analyzed in section 5.1, and the corresponding finite temperature 
parts to be computed in the following. The radiative correction to the renormalized x 
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propagator, Bt{p, x)? has the expression 



Bt{p, x) - Bt=o{p, x) 



1 



+ 



1 1 1 



k^ + X {k + py + xl (exp(-i/?A:o) - 1 



(62) 



where k = (— A^o, k). The ko integral may be performed by making a contour integration 
in the complex ko plane, leading to 



BriP.x) = 5t=o(p,x) + 2 I 



d^k 
(2^ 



k'^ + X 



exp f3^k^ + x]-l 



X 



{k + pj 


^ -P + pI 


{[ 


(k + p) 


-P+pI 


+ Apl [P + x) 1 



(63) 



As we have shown above, in order to study the minimization condition up to next 
to leading order in it is sufficient to study the behaviour of for % — 0. From 
the expression of B^ip, x), doing an expansion in the neighborhood of x = at finite 
external momentum, we obtain 



Bt{po,P, x) = BTipo,P, X = 0) + Kj^ + 0{x), 



(64) 



where i^' is a coefficient of order one, which can be evaluated by analyzing the infrared 
divergences associated with OBt/Ox as x 0. From Eq.(|63D it is possible to prove that 

^^t(p, x) ^ 1 (..^ 

2n/3ipl + f) ^ ' 

as X and, hence, K = — l/27r. Furthermore, the last term in the Eg. (|6^) involves 
higher order contributions in x which are not relevant to solve the minimization condi- 
tion. This is due to the fact that, from the gap equation up to leading order in 1/N, 
Eq.(ig), 

'^^^'^ '^^^ -A.PVX. (66) 



47r/5Vx 



902 1 + 2Att/3^/Xt 
Therefore, the derivative operator involved in the minimization condition, Eq. (^Tj) , only 
receives contributions from the infrared dominant part of dVn.i.o/ dx- 
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An analogous procedure may be used to analyze the behaviour of the vacuum polar- 

(67) 



ization contribution Iljf. 



+ 



-00+.. 27r J (277)3 y^k^.^xm + pf + x] 
{2k + pY{2k + py 



-4 



k^ + X (exp(-2/?A;o) - 1) 



k"^ + X {k + pY + X 

At finite temperature it is easy to demonstrate the transversality of Ui^ by making use 
of the po quantization, pq = ujn = 27m/ jSo- However, general covariance is explicitly lost 



and, hence, the vacuum polarization takes the general formpO 



nr(p,x)=nf(p,x)U'^^ 



p2 



+ n^(p, x) s'^ 



{61 



Recalling the general form of the gauge boson propagator. Eg. (|5^) , it is straightforward 
to show that 

1 2 

el {U^{p,x) + Il^{p,x)] 
1 



det(-zD;J-egnJ, 



1 



[Drip', 
1 



-e^nf(p,x) 



Drip' A) 



^L(p^x, 



(69) 



where, from Eq.([T3), we have 



1 



Dr{k^<P) 
1 

Diik^xJ) 



e + X) + ael^'x 



Therefore, 



det (-iD-l - e^njj 



a (fc2 + x) 



e + el^'-el (nf(fc,x) + n^(A;,x) 
e + el^'-ellil{k^x)\ 7^. 



(70) 



(71) 



This expression is equivalent to that one found by Carrington in the limit of small 
external momentum (po = 0,p ~ 0), while computing the ring diagrams contribution 
to the effective potential in the Landau gauge. Observe that, at zero temperature. 



T=0 



n 



while n 



T=0 



and all logarithmic divergences may be absorbed in the 
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definition of tlie renormalized gauge coupling, wliicli implies that U^^^ = nT'=o, turns 
into its renormalized expression nf"° = TIt=o- Consequently, defining as the finite 
temperature renormalized quantity, the logarithm of the determinant takes the form 



lndetHD-J(fc,0,x) - eluUk,x)]=2\n 



+ In fc^ + e^^-e^nf 



+ In 



(72) 



k"^ + X 

The first two terms in the right hand side of the equation above are gauge independent 
and, for Ilf = 0, they give the well known one loop contribution to the effective potential. 
The last term contains all the gauge dependence and its contribution to the effective 
potential only vanishes in the Landau gauge, a = 0. We shall return to the issue of the 
gauge dependence in section 5.3. 

Since 11^^,, Eq. (p8|) , has only two independent components, we can obtain the func- 
tions nf by computing the expressions of IIoo and 11^ . In fact. 



n 



00 



n 



T=0 
00 



3 (n[ - n[=°) + 2n^. 



(73) 



Moreover, using Eqs. (|49|) and (|67D, we have 



n°?(fc,x)-n°io(A;,x) = h{k,x)-2h{k,x) (74) 
n^'^(fc,x)-n^to(fc,x) = -^l2ik,x)-{k' + ^x){BTik,xJ)-BT=oik,xA)), 



with 



(2fco + Po) 



^+*^ dko r d^k 

<^+ie 27T J (27r)3 1 (^2 + ^) Wj. + + ^ 



+ 



{2ko - poY 



k^ + X 



and 



h{p, X) 



°°+«^ dk 



oo+ie 



k + p) +X 



d^k 



1 



{exp{—i(3ko 



1 



(75) 



2n J (27r)3 (P + x) {exp{-i(3ko) - 1) ' 



(76) 
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The integrals Jj may be computed by performing a contour integration in the complex 
ko plane. We obtain, 

. ^ . d^k 1 
h 



X 



exp l3^Jk^ + x]-l 





P + x)] 


pI+ {k + pf -P 


+ 8pl {P + x) 




Po+ (k+p) - P 


2 

+ Apl (P + 



(77) 



and 



d^k 



1 



(27r)= 



/c2 + X exp /3V + X - 1 



(78) 



It is first interesting to analyze the behaviour of Uipo, p) at x = 0, po = and \p\ <^ T. 
The integral I2 has already been considered in section 4, and its high temperature 
expansion result, Eq.(^), in the limit we are studying, gives /2(x = 0) = 1/12/3^. 
Furthermore, from Eq. ([77|) , it is straightforward to show that 

1 



/i(x = 0,po = 0,p^O) = -- / d\k\- 

The above results imply that, in this momentum regime, 

1 



1 



(79) 



T=0 



-n 



3(3' 



Observe that this contribution gives a static mass term to the abelian gauge field. In 
fact, recalling Eq.(|7lD we obtain that, since nf and come with opposite signs, only 
the longitudinal mode receives a contribution to the mass term equal to 

^2 



Ml 



(81) 



which is the well known Debye screening mass term ^0|- has been recently pointed 
out in the literature that this term is important in order to determine the 

strength of the phase transition within the framework of the improved one loop effective 
potential. Within such approximation, the phase transition appears to be first order 
but, due to the Debye screening suppression of the longitudinal mode, its strength is 
about two thirds of the one obtained using the naive one loop approach. We shall now 
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demonstrate that, when using the self consistent 1/N expansion up to next to leading 
order in 1/N, and computing the value of Ui^ at finite external momentum, extra terms 
appear in the effective potential, which change not only the strength but also the order 
of the phase transition. 

5.3 Minimization of the Effective Potential 



The full effective potential up to next to leading order in 1/N, Eq.(p^ may be rewritten 



as 



Vtree + VJ.^. + - / In 
2, Jk 



Z Jk ^ '- 

+ In + 6^02 _ g2nT 



\e + x) (l + A^Tfefe)/3)+A0V3 
k"^ + X 

- (nf + n^) 



(82) 



where e and A are the renormalized gauge and quartic couplings, respectively, and 
and nj" are the full, temperature dependent radiative corrections contributions with their 
logarithmic divergences removed through the definition of the renormalized couplings, 
in the way explained above. Moreover, we can rewrite the gauge dependent term, V^.d., 
coming from Eq.(^) like 



Vg,a. = \j^ [in (k^ + i??) + In (k^ + i?^) _ (k^ + ^ 



where 



(83) 



(84) 



Similarly to the computation of Bt-, from Eq.(^7D it follows that, in order to evaluate 
the minimum condition up to next to leading order in 1/A^, we only need to study 
the behaviour of TlJ for x — ^- As we have just remarked, the computation of the 
vacuum polarization expression for finite external momentum is crucial in obtaining the 
correct behaviour of the effective potential at small values of x- Recalling Eqs.(|^),(p7|), 
(|7|),(|73) and (^, it is straightforward to find that, for x ^ 0, p'' 7^ 0, 



n^^(po,P,x) = n^^(po,P,X = 0) + 



2ti(3 



+ h.o.{x), 



23 



and 



n!f"(po,P,x) = n!*(po,P,x = o) + 



p2 



+ h.o.ix)- 



86) 



2nP ipl+iP) 

Note that, as implied by transversality, the expression of H^q vanishes in the hmit p ^ 
for Pq 7^ 0. From the above expressions it follows that, at finite momentum transfer and 
in the small x limit, the vacuum polarization components read. 



nf(po,p,x) 
n^(po,p,x) 



nr(po,P,x = o) + ^ + /i.o.(x) 

n^(po,p,x = 0) + h.o.{x)- 



(87) 



It is important to emphasize that, this is the correct behaviour for x ^ at any 



finite \p\ = ypg + which is the relevant regime for the computation of the integrals 
contributing to the minimum of the effective potential. Using the above expressions, 
Eqs.(|6^ ) ,(|87| ) , we can rewrite the minimization condition as 

d dx d 



X 



+ 



+ 



+ 



1 

'n 

In 



'1^ 



In 



502 dx, 

+ 6^02 _ e2 
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g.d. 



nj{k,x = o) + 



vr 

2n8 



e + e^02 - e' (fc, X = 0) + (fc, X = 0) + 

(fc2 + x) (3/A + Brik, X = 0) - ^/{27:^e)) + 

A;2 + X 




The first term in Eq.(p8D involves all the possible gauge dependent contributions to the 
minimum of the effective potential. However, explicitly applying the total derivative to 
the logarithmic expression of K/.d., 



d dx d 
+ - 



dx, 









) Vg.d. 


.^0 2I 


Ikk^ _ 



d 



+ 



dx d 



d(j)^ 

(rI + RI-X 



9) 



502 dx 

and using the fact that {R\ + R^ — x) = 0, we observe that Vg,d. does not give any 
contribution to the minimization condition up to this order in ( Observe that, since 
we are working with regularized integrals with a finite cutoff, we can safely interchange 
derivatives by integrals in the expression above.) Therefore, although the expression 
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of the effective potential is gauge dependent, the minimization condition gives a gauge 
independent result. 

Furthermore, from Eg. (|88D , it comes to notice that, as we discussed above, for the 
minimization of the effective potential, the only relevant contributions from the finite 
temperature radiative corrections come from the infrared dominant terms of Bt{p,x) 
and tifip, x)- III addition, recalling the fact that dxldcj? yx47r/? as x — ^ 0, it readily 
follows that the minimization condition reads, 



X 



+ 



+ 




1 - '^^d{^)/dx 



F + e202 _ e2 (^nf (A;, x = 0) + n^(A;, x = 0) 
1 - '^Vxd{Vx)ldx 



k2 + e202 _ e2 (^nf (A;, x 
l-2^d{^)/dx 



(3/A + 5t(A;,x = 0)) + 




(90) 



The above equation implies that, at next to leading order in the only minimum is 

at 

X™n(0) = 0. (91) 

Thus, the value of xi^P) the minimum remains unchanged after the inclusion of the next 
to leading order contributions. As we discuss in section 4, Eg. (pT]) leads to a nontrivial 
solution for the vacuum expectation value of the scalar field, = — 1/(6/3^), for 
temperatures below the critical temperature, = 6v'^. Hence, the vacuum expectation 
value of the scalar field varies continuously at the critical temperature, characterizing a 
second order phase transition. 

6 Comparison with the One Loop Approximation 

The results of section 5 show that, in the presence of a large number of scalars and in 
the region of parameters e^/X <^ N, in which the scalar loop contributions are enhanced 
in comparison to the gauge loop ones, the phase transition is second order. In spite 
of that, the gauge contributions to the effective potential at next to leading order in 
resemble those ones contributing to the effective potential in the improved one 
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loop approximation. In fact, for small (f), and at T = Tq, at which the quadratic term 
in (j) vanishes, one would expect that due to the presence of the gauge contributions, 
a dominant negative cubic term would appear, leading to a nontrivial minimum in the 
effective potential and, hence, to a first order phase transition. However, the fundamental 
difference comes from considering the finite external momentum contributions to the 
radiative corrections at finite temperature. 

Quite generally, as we have already said, the integral 

It = 1 f ln(A;2 + m') (92) 

2 Jk 

may be evaluated in the high temperature regime to give, 

'^-'r^.^^-^^^k.o.i'nV). (93) 

Therefore, since the zero temperature gauge boson mass is equal to — e^cp one would 
naively expect that the inclusion of the gauge field would induce the generation of a 
quadratic term in 4>, which would modify the definition of the transition temperature, as 
well as the appearance of a cubic term, which would lead to a first order phase transition. 
In the next to leading order in 1/N analysis, these two effects are not present, and, 
hence, some sort of finite temperature screening, which modifies the nature of the phase 
transition, should occur. One type of screening, which is already contained within the 
improved one loop approximation, is that one induced by the generation of an effective 
temperature dependent mass for the longitudinal gauge boson mode. As we already 
discussed above, when compared with the naive one loop result, the Debye screening 
causes the reduction in two thirds of the value of the coefficient of the cubic term. In 
our analysis, however, the preservation of the value of the critical temperature and of 
the nature of the phase transition after considering the next to leading order in 1/N 
effects is a result that would be obtained even in the case of a cancellation of the 
independent radiative corrections contributions to the photon and scalar propagators. 
In fact, within the 1/N expansion, there is a different source of screening, which can 
only be obtained by considering the dependent, finite momentum contributions to the 
radiative corrections. 
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For values of (p close to the minimum of the effective potential (x — 0), and at 
momentum ^ x, the evaluation of the finite temperature radiative corrections to 
the photon and scalar propagators have as outcome the replacement of the functional 
dependence on cfP' by a dependence on the combination [cf? — ^Jxl (27r/5)]. In fact, as it is 
clearly seen from Eqs.(p2D, (|85| ) -(|87| ) , in this momentum regime, the finite temperature 
effective mass of the transverse gauge boson mode is given by 



r 



(02-v^/(27r/5))+/i.o.(x) 



(94) 



To easily understand the relevance of this effect, let us first analyze the behaviour of x 



at T = Tr 



2A-k(3, 



+ 



( At. 



+ 



6 



(95) 



Expanding Eq.(^) in the neighborhood of the origin we obtain 

, \ 2' 



At„ 



(96) 



From Eq. (|96|) it follows that at T = Tc and for the momentum regime under considera- 
tion, the effective mass of the photon transverse mode in the neighborhood of the origin 
is given by 



r 



A 



(97) 



Thus, since for the same range of parameters. 



X 



(98) 



the effective squared mass of the photon transverse mode behaves like rather 
than the expected 0^ dependence. Within the region of parameters considered above, 
the radiatively corrected squared mass of the a particle also behaves as 

The dependence of in the neighborhood of the origin and at T = Tc has two 
effects. First, it explains the cancellation of the quadratic term in at T = and 
hence the preservation of the value of the critical temperature at next to leading order 
in Second, it implies that, in the neighborhood of the origin, and independently 
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of the existence of - independent radiative corrections contributions, no cubic term is 
induced and thus, at T = Tc the finite temperature effective potential is well described 
by a positive, quartic term. 

It is important to remark that the expressions we obtain for the radiatively corrected 
photon and scalar propagators are only valid for T, k"^ ^ x- Moreover, the cubic term 
arising in the integral. Eg. (p2D , receives the most relevant contributions from the mo- 
mentum regime = 0, k^ = 0{m'^). Thus, the present analysis would fail if in the 
neighborhood of the origin the relation 0^ ^ x were not fulfilled. However, as readily 
seen from Eg . (|98|) , x ^ 0^ and the effective potential close to the origin is well described 
by our approach. 

7 Conclusions 

In this article we have analyzed the finite temperature phase transition of scalar elec- 
trodynamics with N flavors of complex scalar fields, by means of a large N expansion, 
up to next to leading order in We have shown that the effective potential takes 

a compact and simple form, which, although gauge dependent, leads to gauge invariant 
results for the critical temperature as well as for the extrema of the potential, which are 
well defined physical guantities. At leading order in the effective potential coin- 

cides with the one of the 0(2N) vector model, already studied in the literature. At the 
critical temperature = 6f^, the system develops a symmetry restoration phase tran- 
sition, with an order parameter which varies continuously from zero to nonzero values. 
The phase transition is hence second order. 

At next to leading order in the expansion, the effective potential receives con- 
tributions from the radiatively corrected a and gauge field propagators. The gauge field 
contribution is essentially that one found at one loop, but with the photon corrected by 
vacuum polarization effects. The (p dependence of the vacuum polarization effects is cru- 
cial in order to define a correct expansion. We have shown that, when the radiative 
corrections are considered at finite external momentum, their effects lead to an effective 
screening of the gauge boson mass, which change the nature of the gauge boson loop 
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contributions to the effective potential. Due to this screening effect, no extra 0^ term is 
induced in the neighborhood of the origin at the leading order critical temperature and, 
therefore, the critical temperature remains the same as in the leading order analysis. 
More generally, the position of the minima of the effective potential is not modified by 
next to leading order effects and, in addition, no new minimum of the effective potential 
appears at next to leading order. Hence, the order parameter has the same behaviour 
as in leading order in 1/N and the phase transition remains second order. 

The results of this work show that in the presence of a large number of flavors of 
complex scalar fields, the symmetry restoration phase transition of scalar QED becomes 
second order. Within the context of this study, however, we can not rule out the presence 
of a critical value of A^, below which the transition is first order. In general, we expect 
the range of validity of the large N approximation to depend on the relation between the 
gauge and quartic couplings. More specifically, the expansion considered in the present 
work assumes the dominance of the scalar loops contributions to the effective potential. 
Therefore, it is rigorously vahd for values oi e^/X <^ N. Thus, for low values of N, the 
expansion becomes more rehable for a scalar Higgs mass, = A0^/3 larger than the 
gauge boson mass, rri^ = e^0^, ml/rri^ > 1, for which the fluctuations of the scalar fields 
are enhanced. 
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FIGURE CAPTIONS 



Fig.l. Multiloop diagrams which contribute to the effective potential at next to leading 
order in Sohd and dashed hues denote the propagators of the scalar fields, 0a,i and 
X, respectively, while the curved lines denote the gauge field propagator. 

Fig. 2. Same as in Fig. 1, but considering those diagrams which involve a mixing of 
X and Af,. 
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